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Abstract 



Radiation (CBR) anisotropy depends upon the present cosmological param- 
eters as well as inflationary parameters. The analysis of a sampling- variance- 



^ , limited CBR map offers the most promising means of detecting tensor per- 

H ; 

5t , turbations, but will require evaluation of the predicted multipole spectrum 

for a very large number of cosmological parameter sets. We present accurate 
polynomial formulae for computing the predicted variance of the multipole 
moments in terms of the cosmological parameters Q/^, HqH"^, VLsh"^-, N^, and 
the power- law index nx which are accurate to about 1% for / < 50 and to 
better than 3% for 50 < / < 100 (as compared to the numerical results of a 
Boltzmann code). 

PACS index numbers: 04.30.+X, 98.80.Cq, 98.80.Es 



I. Introduction 



Inflation makes three robust predictions: flat Universe and nearly scale-invariant spectra 
of scalar (density) and tensor (gravitational-wave) metric perturbations The scalar 
and tensor |Q perturbations arise from quantum mechanical fluctuations on very 
small scales during inflation and are stretched to astrophysically interesting scales by the 
tremendous growth of the cosmic scale factor during inflation. 

Both the scalar and tensor perturbations give rise to anisotropy in the temperature 
of the Cosmic Background Radiation (CBR) seen on the sky today, most conveniently 
described by their contribution to the multipole decomposition of the CBR temperature 

oo / 

STie,(P)/T = J2 E aimYUOA)- (1-1) 

1=2 m=-l 

The scalar and tensor contributions to the anisotropy predicted by inflation are uncor- 
related and statistical in character. The individual multipoles that describe our sky are 
given by the sum of a scalar plus tensor contribution, with these contributions being 
drawn from gaussian distributions with variances and which are related 

to the properties of the inflationary potential and cosmological parameters. Because the 
scalar and tensor contributions are uncorrected, the variance (|a/mP) = (khnP) + (kimP)- 
The expected scalar 0] and tensor contributions are shown in Fig. 1 for a nominal set of 
cosmological parameters. 

CBR anisotropy offers a very promising means of testing inflation as well as deter- 
mining the scalar and tensor perturbations. If both can be measured, then information 
about the underlying inflationary potential can be derived (value of the potential and its 
first few derivatives at a point) ||^. Key to doing this is a high-angular-resolution (better 
than 0.5°), sampling- variance-limited map of the CBR sky. (Because there are only 21 + 1 
multipoles, sampling variance limits the accuracy to which (|a/mP) can be measured - a 
relative precision of ^2/(2Z + 1).) Three proposals have been made to NASA (FIRE, PSI, 
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and MAP) for a satellite-borne experiment and another to ESA (COBRAS/SAMBA). A 
satellite could be launched as early as 1999. 

The separation of the scalar and tensor contributions to CBR anisotropy is likely 
to be done by maximum likelihood techniques and will require accurate predictions for 
the scalar and tensor contributions (1% or better) to the anisotropy for many sets of 
cosmological and inflationary parameters. At present, achieving such precision requires 
numerically integrating Boltzmann equations, which is very time consuming (typically 
requiring many hours on a powerful workstation for one set of cosmological parameters) 
0. The need for a fast and accurate approximation scheme is manifest. 

Many analytic approximations to the tensor angular power spectrum have been ex- 
plored [0 . The most accurate is as time consuming as numerically integrating the Boltz- 
mann equations 0, and even schemes with less accuracy require significant computation 
time (tens of minutes). This motivated the present approach - a polynomial fit that can 
be evaluated very rapidly (much less than a second for a set of cosmological parameters). 

The tensor angular power spectrum, Ci{P) = depends upon a set of cosmo- 

logical and inflationary parameters denoted here by P. The set P is: the baryon density, 
fis/i^; the matter density, VLqK^] the level of radiation in the Universe, parameterized by 
the equivalent number of massless neutrino species A^;^; the present vacuum-energy den- 
sity, = 1 — ^o! Ei'iid the primordial power-law index of the tensor perturbations, 
riT {fiT = for scale invariant tensor perturbations). We will only be concerned with 
the shape of the angular power spectrum; the overall amplitude of the angular power 
spectrum, conveniently specified by C2, depends upon the inflationary parameters (the 
value of the inflationary potential in Planck units) as described elsewhere 0. In most 
approaches, the shape {C1/C2) and the overall amplitude (C2) are determined indepen- 
dently. 

The dependence upon the parameters is simple to explain. The shape of the angu- 
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lar power spectrum depends on the redshift of last scattering and the evolution of the 
gravitational waves after they enter the horizon, which depends upon the evolution of 
the cosmic scale factor. The redshift of last scattering (more precisely, the peak of the 
visibility function; see Appendix B) depends upon the baryon density, the matter den- 
sity, and the level of radiation in the Universe. The evolution of the scale factor of the 
Universe around last scattering depends upon the relative levels of matter and radiation 
through the value of the scale factor at matter-radiation equality. 



The recent evolution of the scale factor, which depends upon as well, is also important 
for gravitational wave modes which have recently entered the horizon and influence the 
low / multipoles. Finally, the primordial spectral shape of tensor perturbations (described 
by tit) also affects the shape of the angular power spectrum. 

We have engineered our fits based upon the numerical results of the Boltzmann code 
written by Dodelson and Knox [|ry]; they believe that their code is accurate to better than 
1%. We expand Ci{P) around a fiducial set of parameters: for f^A = {VIq = 1), h = 0.5, 
A^^ = 3, Qeh^ = 0.0125, and ut = 0; and for 7^ 0, Qoh^ = 0.125, Qeh^ = 0.0125, 
^A = 0.5, Ni, = 3 and tit = 0. These two cases are treated in the next two Sections. We 
end with a brief discussion of the accuracy of our fits. 



II. Zero Cosmological Constant 

For ^A = 0, the cosmological parameters are [h, N^,QBh'^ j^t)- We define a parameter 
vector, P = {[hf{N^)],N^, [nBh'^]-\nT), where 




(1.2) 




(2.1) 



which is related to Req of Eq. ( ]1.2| ) . 
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We expand the tensor multipole spectrum for P around the tensor multipole spec- 
trum for P = Q = (0.5, 3, 80, 0) [which corresponds to h = 0.5, = 3, Qsh^ = 0.0125, 
and riT — 0]- We write 



— [l + g{L)nT\ 



C2{P) C^iQ) V2, 



l + Eti/i(^i-a,2) 



(2.2) 



where 



/i(Pi-gi,0 = E«i(0 (Pi-QiY 

f2iP2-Q2,l)^j2bjil) {P2-Q2y 

h{Pz-Qz.i) = jZdj{i) (Ps-QsY 



(2.3) 



The coefficients are found numerically by fitting to the variation due to each parameter 
separately. Identifying the relevant vector of parameters minimizes the need for cross 
terms. 

The Pi coefficients are 



ai(/) 

02(0 
asil) 



-0.4025 (—\ - 0.3375 ( —\ - 2.3441 ( + 2.0125 ( —\ 
\100 [lOO \lOO \100 



1.1271 ( — I - 2.0614 ( — 

UOO/ uoo 



I + 10.9105 (-'—] - 10.8102 (—\ 
' I 100/ UOO/ 



-2.7875 (-'—] + 10.9417 ( - 36.4181 ( —\ + 41.1617 ( 

UOO/ UOO/ UOO/ UOO/ 



aM) = 3.4719 f—^ - 16.8888 f—^ + 51.2793 f—^ - 61.1775 f—^ 
UOO/ UOO/ UOO/ UOO/ 



(2.4) 



The P2 coefficients are 

bi{l) = 10-3 



-0.0361 I 1 +5.3684 ^ 



100, 



, + 3.3680 

100/ UOO, 



-2.0872 (—^ + 1.2639 (— 
UOO / I 100 
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HI) 



10" 



-0.1105 



100, 



+ 0.1237 



100, 



1.8853 



100, 



+2.2411 



100 



- 1.0243 



100 



10- 



0.1297 



100, 



0.5125 



100^ 

\ 4 



+ 2.3069 — 

uoo. 



(2.5) 



The A coefficients are 



di{l) = 10"^ 
d2{l) = 10-^ 

^3(0 = 10-^° 



-0.1972 



100, 



12.3357 



100, 



- 0.7469 (—\ + 4.0362 (— 
llOO / 1 100 



-0.0241 (—\ - 1.7123 (—\ - 1.1385 ( + 0.6791 (— 
llOO/ I 100/ llOO/ I 100 



0.2028 ( 1 + 12.3878 I — 
I 100 / V 100 



I + 8.9878 f-^^ - 6.1211 f-^^ 
' I 100 / V 100 / 



(2.6) 



The coefficient involving ut is 



^(0 



-0.4764 



1 _ g-{«-2)/30' 



+ 1.6734 



1 _ g-2(/-2)/30 



- 4.0400 



+4.6345 



1-e 



-4(i-2)/30 



2.2942 



1-e 



-5(;-2)/30 



1 _ g-3(«-2)/30 

(2.7) 



III. Nonzero Cosmological Constant 



For Q\ > 0, the cosmological parameters are (ilo^^, ^sh^, ^a, -^i/, ^t-t)- The procedure 
of fitting is similar to the Qa = case, but shghtly more comphcated because one "cross 
term" is required to achieve sufficient accuracy. 

Here we define the parameter vector, P — f{Ni,)^/n^ , A^,^, Qa, ^t-t)- 

Again, we expand the C^'s for P around the C;'s for a standard set of cosmological 



parameters, S. For Qa < 0.36, we take S — (0.5, 3, 80, 0, 0) [which corresponds to 
h = 0.5, = 3,nBh^ = 0.0125, i^A = 0, and ht = 0]. For 0.36 < I^a < 0.8, wc take 
S = (0.3536, 3, 80, 0.5, 0) [which corresponds to h ^ 0.5, = 3, O^/i^ ^ 0.0125, Oa = 
0.5, and ut — 0]. As before, we write 

Q{p) _ci{s) (iy\,^ 1 / l + T.Uf^{p^-s^^)+v{lm-Slm-s,) 



C2{P) C2{S) 12, 



1 + Eti MP^ - S^, 2) + r]{2){P, - S,){P, - S,) 



where 

^(Pl-S^,l) 

/2(P2-^2,0 

fi{Pi — Si, I) 
For Qa < 0.36, 



^A,(0 (Pi-^i^ 
X:i^, (0 (Ps-^ay 

Ei?,(0 (^4-^4^, 



A,(0 = %(0, 5,(0 = 6,(0, i^,(0 = c?,(0; 

K0 = 5(0, ^(0 = 0. 



The P4 coefficients are 



£;i (2) = -6.0163 X 10-^ £;i(3) 



^2(2) = -3.2853 X 10"^ ^2(3) 



E^{2) = -1.8426 X 10-^ E^{?,) 
^4(2) = -2.8162 X 10"^ ^4(3) 



-1.4845 X 10-^ 
-8.4852 X 10"^ 
-4.0997 X 10-^ 
-8.6348 X 10-^ 



Ei{l>A) = 0.1 



(3.2) 



(3.3) 



£^1(4) 
^2(4) 
£^3(4) 
^4(4) 



-4.3594 X 10~^ 
-2.5631 X 10"^ 
-9.4267 X 10"^ 
-3.1561 X 10"^ 



0.1483 i-'—] - 2.4098 { —\ + 0.3211 {— 
1 100 y I 100 / I 100 



- 0.6630 



100 



E2{1 > 4) 
Es{l > 4) 
E4(l > 4) 



0.1 



0.1 



0.1 



0.0768 



0.0580 



0.0560 



/ 



100^ 

100^ 

Too 



1.5272 ( — 1 - 0.2895 ( — ) + 0.1754 ( — 
1 100/ 1 100/ 1 100 



- 1.0585 ( — I - 0.1897 f — ) + 0.2662 ( — 

Vloo; viooy vioo 

- 1.6921 ( —] - 0.8962 ( —\ + 1.0897 (— 

I 100 / 1 100/ I 100 



All other coefficients have been given in the previous Section. 
For 0.36 < Q\ < 0.8, the Pi coefficients are 



Ai{l) = -1.0313 



I \ „„_ / / \^ f I V 



100 



1.0270 



100 



7.6262 



100 



+6.8176 (—] - 2.9419 ( 

Viooy Vloo; 

A2(l) = 3.5896 f—^ - 11.9704 f—^ + 37.1339 f—^ 
I 100/ 1 100/ llOO/ 



Ml) 



-7.6185 



100, 



-35.5609 



+ 34.4566 



+95.4331 



—] + 19.0353 f-^^ 
100/ llOO/ 



/ 



100 



92.1698 



/ 



100, 



—] - 53.5697 f-^ 
100/ llOO 



/ 



/ 



^4(0 = 6.5272 — - 32.9224 + 84.729^ , 

I 100/ I 100/ VIOO 



-91.4209 ( —\ + 51.8163 (— 
llOO/ I 100 



The P2 coefficients are 

= 10-2 



0.1629 (-^ 1 + 1.4392 I — 
V 100 / V 100 , 



- 2.2268 



100, 



^2(0 = 10 



-4 



+2.8450 f-^^ - 1.5396 f-^ 
1 100 / 1 100 , 



0.3608 (-!—] - 6.3470 ( —\ + 3.4559 (— 
\lOO I 100/ I 100 



(3.5) 
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10" 



0.1672 ( I - 0.8800 f — | + 2.9272 f — 

1 100/ liooy llOO 



-3.6758 (—] + 1.6453 ( — 
I 100/ I 100, 



(3.6) 



The Ps coefficients are 



DiiD 

D2{1) 



10" 



-0.1039 (—\ + 9.2599 ( 4^-] + 0-1085 (— 
I 100 / 1 100/ I 100, 



+ 2.4910 



100 



10 



10" 



-7 



-0.2716 l-i-] - 16.4608 ^ 



100 



I - 11.4085 (—^ + 8.5966 (— 

100/ I 100/ llOO 



0.0962 



100 



+ 2.8879 



100, 



+ 3.0975 



100 



- 3.0516 



100 



(3.7) 



Tlie P4 coefficients are 



E,i2) = 
E2(2) = 
^3(2) = 
E^{2) = 

Ei{l > 4) 
E2{1 > 4) 
Es{l > 4) 
E^(l > 4) 



-1.3087 X 10~\ 
-1.1277 X 10"\ 
-1.3674 X 10-\ 
-6.8693 X W-\ 



El (3) 
E2(S) 
^3(3) 
£^4(3) 



-3.1368 X 10-2, 
-2.7633 X 10-2, 
-3.4105 X 10-2, 
-1.7588 X 10-^ 



8.4146 X 10"^ 
7.3835 X 10-^ 
8.9756 X 10"^ 



El (4) = 

^2(4) . 

^3(4) = 

^4(4) = -4.5389 X 10 



-2 



0.1 



0.1 



0.1 



0.5211 

0.4547 

0.5924 
/ 



/ 



100^ 
100, 



- 5.7831 



- 5.2202 



- 7.3578 



4) + 1.0999 



/ 



,100. 
,100, 

V 2 



- 0.1399 

- 1.5359 



,100, 

.100, 

,100, 
\ 3 



- 1.2519 

- 0.3364 

+ 1.3647 

/ 



,100, 
,100, 
,100, 

, 4 



0.5016 - 7.0197 - 0.2120 +2.9141, 

I 100/ llOO/ llOO/ VlOO 



(3.^ 
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The coefficient involving ut is 



-0.4670 



1 _ g-a-2)/30' 



+ 1.7234 



1 _ g-2(«-2)/30 



4.2724 



-4.9757 



1 - e 



-4(«-2)/30 



- 2.4571 



1 - e 



-5(/~2)/30' 



1 _ g-3(Z-2)/30 
(3.9) 



Note that u{l) ~ g{l). 

The cross-term coefficient is 



^(0 



-3.7201 f— ^ + 29.9770 f— ^ - 117.6372 f — ^ + 236.1950 f— ^ 
I 100/ I 100/ I 100/ I 100/ 



-257.9630 f — ^ + 147.3101 f— ^ 
V 100 / V 100 / 



- 34.5437 



100 



(3.10) 



IV. Discussion 



By identifying the relevant cosmological parameters we have developed a fast (<^ 1 sec) 
and accurate (few percent or better) algorithm for computing the shape {C1/C2) of the 
tensor angular power spectrum for a primordial tensor power spectrum of the form 



PT{k) oc {kroY^k-^, 



(4.1) 



where tq is the conformal time today. Our algorithm employs a polynomial in the pa- 
rameters f^A, ^oh"^-, ^sh"^, TT'Ti and Nij. 

To assess the accuracy of our algorithm we sampled the following parameter intervals 
uniformly and at random: 0.35 < h < 0.8, 2 < < 12, 0.005 < Vteh^ < 0.03, 
—0.3 < "H-r < 0, and < i^A < 0.8. In Figs. 2 and 3 we show the histograms of 
the maximum error in C1/C2 for ^a = and 7^ respectively. For i^A = the 
maximum error (relative to the Boltzmann code of Ref. [0) is less than about 0.5% for 
/ < 50 and less than about 2% for / < 100. For i^A 7^ 0, the accuracy is slightly worse, 
better than about 1% for / < 50 and better than about 3% for I < 100. In the case of 
Q\ the largest errors occur for fig/i^ < 0.05 (large f^A and small h). 



9 



The tensor contribution to the quadrupole plays a special role. It provides a con- 
venient overall normalization for the angular power spectrum and can be related to the 
value of the inflationary potential when the comoving scale = Hq crossed outside the 
horizon during inflation M: 



K/m^i = 0.66 [1. - (/i,') + 0.1)nT\ C^/ f^'> , (4.2) 
where the functions fj!''^\^ljy) are given by 

f^\n^) = 1. - o.osr^A - o.i^i 

/W(Q^) = 0.58 - O.SOl^A + O.Sl^i - 0.88(]i. (4.3) 

The dependence of this relationship on cosmological parameters other than is much 
less significant 

There is even more motivation for developing a fast and accurate algorithm for the 
scalar angular power spectrum. However, this task is more challenging: the power spec- 
trum has more structure (cf., Fig. 1) and that structure extends to higher multipoles. 
We are currently working on an algorithm for the scalar angular power spectrum. 
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Appendix A: Table of standard C/'s 



I 


i{i+i)Ci{nA=o) 

6C2(f2A=0) 


i{i+i)Ci{nA=o.5) 

6C2(nA=0.5) 


2 


1.000000000 




1.000000000 




3 


7.85882861E- 


•01 


8.12787313E- 


■01 


4 


7.43475686E- 


•01 


7.75447410E- 


•01 


5 


7.35658928E- 


■01 


7.69590303E- 


■01 


6 


7.38104089E- 


■01 


7.73328636E- 


■01 


7 


7.43851511E- 


■01 


7.80190352E- 


■01 


8 


7.50438650E- 


•01 


7.87827791E- 


■01 


9 


7.56898572E- 


•01 


7.95302374E- 


•01 


10 


7.62844076E- 


■01 


8.02234328E- 


■01 


11 


7.68121692E- 


■01 


8.08471560E- 


■01 


12 


7.72700749E- 


•01 


8.13985567E- 


•01 


13 


7.76567875E- 


•01 


8.18762161E- 


•01 


14 


7.79782672E- 


•01 


8.22865837E- 


•01 


15 


7.82327767E- 


■01 


8.26275290E- 


■01 


16 


7.84314505E- 


■01 


8.29111207E- 


■01 


17 


7.85671940E- 


•01 


8.31292902E- 


•01 


18 


7.86573958E- 


•01 


8.33010933E- 


•01 


19 


7.86864814E- 


•01 


8.34090814E- 


•01 


20 


7.86805032E- 


■01 


8.34819817E- 


■01 


21 


7.86129624E- 


■01 


8.34902088E- 


■01 


22 


7.85208901E- 


•01 


8.34747737E- 


•01 


23 


7.83648 11 5E- 


•01 


8.33916120E- 


•01 


24 


7.81949251E- 


•01 


8.32965076E- 


•01 


25 


7.79570004E- 


■01 


8.31288581E- 


■01 


26 


7.77161094E- 


■01 


8.29611689E- 


■01 


27 


7.74020055E- 


•01 


8.27148151E- 


•01 


28 


7.70959238E- 


•01 


8.24805533E- 


•01 


29 


7.67105366E- 


•01 


8.21604444E- 


•01 


30 


7.63442668E- 


■01 


8.18647318E- 


■01 


31 


7.58920711E- 


■01 


8.14753604E- 


■01 


32 


7.54699914E- 


•01 


8.11226404E- 


•01 


33 


7.49551094E- 


•01 


8.06680893E- 


•01 


34 


7.44812364E- 


•01 


8.02623603E- 


•01 


35 


7.39076276E- 


■01 


7.97464933E- 


■01 


36 


7.33855970E- 


■01 


7.92913085E- 


■01 


37 


7.27571031E- 


■01 


7.87178514E- 


•01 


38 


7.21903458E- 


•01 


7.82165288E- 


■01 


39 


7.15107322E- 


•01 


7.75891111E- 


•01 


40 


7.09025177E- 


•01 


7.70447667E- 


•01 
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I 


fjyi-r -i- t-A — 






6C2(J^A=0) 




6C2(nA=0.5) 


41 


7.01754819E- 


■01 


7.63669098E-01 


42 


6.95290317E- 


■01 


7.57825769E-01 


43 


6.87581696E- 


•01 


7.50576759E-01 


44 


6.80766517E- 


•01 


7.44363318E-01 


45 


6.72655788E- 


•01 


7.36677916E-01 


46 


6.65520307E- 


■01 


7.30122320E-01 


47 


6.57043916E- 


■01 


7.22034754E-01 


48 


6.49618128E- 


•01 


7.15164612E-01 


49 


6.40812380E- 


•01 


7.06709069E-01 


50 


6.33125395E- 


•01 


6.99550707E-01 


51 


6.24026867E- 


■01 


6.90761689E-01 


52 


6.16107568E- 


■01 


6.83341454E-01 


53 


6.06751841E- 


•01 


6.74252272E-01 


54 


5.98629552E- 


•01 


6.66596862E-01 


55 


5.89051404E- 


•01 


6.57240302E-01 


56 


5.80754839E- 


■01 


6.49375955E-01 


57 


5.70988919E- 


■01 


6.39784406E-01 


58 


5.62546297E- 


•01 


6.31736935E-01 


59 


5.52626894E- 


•01 


6.21943006E-01 


60 


5.44065284E- 


•01 


6.13736759E-01 


61 


5.34026679E- 


•01 


6.03773187E-01 


62 


5.25372061E- 


■01 


5.95432062E-01 


63 


5.15247817E- 


•01 


5.85331082E-01 


64 


5.06526007E- 


•01 


5.76878240E-01 


65 


4.96347985E- 


•01 


5.66670612E-01 


66 


4.87585435E- 


•01 


5.58131161E-01 


67 


4.77383343E- 


■01 


5.47844948E-01 


68 


4.68606069E- 


•01 


5.39242990E-01 


69 


4.58408817E- 


•01 


5.28906521E-01 


70 


4.49640932E- 


■01 


5.20265325E-01 


71 


4.39477440E- 


•01 


5.09906642E-01 


72 


4.30741650E- 


■01 


5.01246850E-01 


73 


4.20639523E- 


•01 


4.90893728E-01 


74 


4.11958163E- 


■01 


4.82237697E-01 


75 


4.01942685E- 


■01 


4.71915135E-01 


76 


3.93337874E- 


•01 


4.63283465E-01 


77 


3.83432932E- 


■01 


4.53015076E-01 


78 


3.74925150E- 


•01 


4.44428914E-01 


79 


3.65153926E- 


•01 


4.34239289E-01 


80 


3.56761877E- 


•01 


4.25716630E-01 
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I 


l(l~\-l)c,(n^—o) 




6C2(r!A=0) 


6C2(r!A=0.5) 


81 


3.47146488E-01 


4.15627684E-01 


82 


3.38887429E-01 


4.07185415E-01 


83 


3.29449095E-01 


3.97220535E-01 


84 


3.21338323E-01 


3.88874768E-01 


85 


3.12097094E-01 


3.79055720E-01 


86 


3.04149866E-01 


3.70821411E-01 


87 


2.95121515E-01 


3.61164479E-01 


88 


2.87354633E-01 


3.53060287E-01 


89 


2.78552720E-01 


3.43581079E-01 


90 


2.70980245E-01 


3.35623379E-01 


91 


2.62418690E-01 


3.26338036E-01 


92 


2.55052066E-01 


3.18539155E-01 


93 


2.46743720E-01 


3.09462129E-01 


94 


2.39593392E-01 


3.01833443E-01 


95 


2.31549532E-01 


2.92978236E-01 


96 


2.24624559E-01 


2.85530583E-01 


97 


2.16855038E-01 


2.76910111E-01 


98 


2.10163325E-01 


2.69653750E-01 


99 


2.02676236E-01 


2.61279304E-01 


100 


1.96225065E-01 


2.54224121E-01 



Appendix B: Last Scattering 

In the early Universe, matter and radiation were in good thermal contact, because of the 
rapid interactions between the photons and electrons. As the temperature dropped below 
0.3eV, electrons combined with protons to form neutral hydrogen ("recombination") at 
a redshift of around 1300. With the disappearance of free electrons, the photon mean 
free path became very large (> H~^) and matter and radiation decoupled at a redshift 
of around 1100. Last scattering is crucial in calculating the CBR anisotropy. ||lT| 

The redshift of last scattering, 2;lsS; is given by the peak of the visibility function 
g{z) = e~'^dr/dz; g{z)dz measures the probability that a given photon suffered its last 
scattering in the redshift interval {z, z+dz). The optical depth (measured from the 
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present back to redshift z) is given by 

t{z) = C ^ ne{z) (7t, 



(6.1) 



where dt is the proper time interval, ne{z) is the electron number density, and ax — 
6.65 X 10~^^cm^ is the Thomson scattering cross-section. The electron number density 
depends on flsh^ and H{t), the Hubble parameter at time t. At the relevant times 
(around recombination and last scattering), H{t) only depends on floh'^ and N^. Hence, 
zlss only depends on flsh^, ^oh^, and Ni,. Numerically, we find 



^Lss = 1104.37 + Az(QB/^^ Qq/i^) + A^(Qo^^) + ^z{N^, ^Bh'^, ^oh'^), 



(6.2) 



where 



Az{nBh^,^oh^) = 0.5285 [(Qfi/i^)-^ - 0.0125"^ 



floh 
0.25 



2 \ 0.31 



Az{noh'^) = 73.21 



-7.022 X 10"^ [(Qfi/i^)-^ - 0.0125"^ 

1 2 



12.06 



Qoh^ - 0.5 



Az{K, neh^, noh'^) = 0.3823(iV^ - 3) 



fi^/^n -0.756 



0.25 



2 \ -0.46 



(6.3) 



^0.0125; \0.25 J 
Our fitting formula is accurate to Az = ±1 for the parameter ranges of 0.1 < Qoh^ < 0.64, 
0.005 < Ofi/i^ < Q 2<N^< 12. 
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Figure Captions 

Fig.l The predicted scalar and tensor contributions. The set of cosmological 
eters is: h = 0.5, = 0.05, A^^ = 3, = 0, ut ^ 0. 

Fig. 2 Histograms of maximum error for I < 50 and I < 100 (Qa = 0). 
Fig. 3 Histograms of maximum error for / < 50 and / < 100 (^a > 0). 
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Log 1 



Figure 1: The predicted scalar and tensor contributions. The set of cosmological param- 
eters is: h = 0.5, = 0.05, K = 3, I^a = 0, tit = 0. 
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Figure 2: Histograms of maximum error for I < 50 and I < 100 (Qa = 0). 
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Figure 3: Histograms of maximum error for I < 50 and I < 100 (Qa > 0). 
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